Spin drag in an ultracold Fermi gas on the verge of a ferromagnetic instability 
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Recent experiments [Jo et a/., Science 325, 1521 (2009)] have presented evidence of ferromagnetic 
correlations in a two-component ultracold Fermi gas with strong repulsive interactions. Motivated 
by these experiments we consider spin drag, i.e., frictional drag due to scattering of particles with 
opposite spin, in such systems. We show that when the ferromagnetic state is approached from 
the normal side, the spin drag relaxation rate is strongly enhanced near the critical point. We also 
determine the temperature dependence of the spin diffusion constant. In a trapped gas the spin 
drag relaxation rate determines the damping of the spin dipole mode, which therefore provides a 
precursor signal of the ferromagnetic phase transition that may be used to experimentally determine 
the proximity to the ferromagnetic phase. 
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Introduction. — Interest in electronic transport ranges 
from everyday applications to fundamental physics. One 
of the most interesting phenomena that spans this entire 
range, is the influence of a thermodynamic phase transi- 
tion on the electrical conductivity. The most direct ex- 
ample is the phase transition from a normal conductor to 
a superconductor characterized by a vanishing resistance. 
The applications of this phenomenon are ubiquitous and 
the basic physics that underlies the transition in super- 
conductors, the Bose-Einstein condensation of fermionic 
pairs, has emerged in research fields from astroparticle 
physics [Tj to cold-atom systems [2j [3] . 

A system in between these two temperature extremes, 
in which analogies of superconductivity have been pre- 
dicted, is that of a two-dimensional (2D) electron-hole 
bilayer [4j [5] . In this case the pairs that condense are 
excitons formed by electrons from one layer with holes in 
the other. The relevant transport probe is in this case 
the Coulomb drag measurement [6 : a current / is driven 
through one layer, known as the "active" layer, causing a 
voltage drop Vb in the other. As the layers are separated 
by an essentially impenetrable tunnel barrier, the volt- 
age drop is predominantly caused by Coulomb scattering, 
and the drag resistivity pr> = Vd/I has the characteristic 
quadratic Fermi-liquid-like low-temperature dependence 
Pd cx T 2 . When the excitons undergo Bose-Einstein con- 
densation, however, the drag resistivity is predicted to 
jump from the relatively small value proportional to T 2 
to a value equal to the ordinary resistivity of the active 
layer [7]. Although conclusive evidence of exciton con- 
densation is still lacking, two experimental groups [8j [9] 
have recently reported the observation of an upturn in 
the drag resistivity as the temperature is lowered. This 
upturn is interpreted as being due to strong pairing fluc- 
tuations that precede exciton condensation [10] and thus 
serves as a precursor signal for the transition, similar to 
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FIG. 1: (Color online) Spin drag relaxation rate l/r s d(T) 
as a function of temperature T, for various values of the in- 
teraction parameter k F a. The Fermi energy is denoted by 
sf — k^T F = h 2 k F /2m. Note that for k F a > tt/2, the spin 
drag relaxation rate shows a distinctive upturn when the crit- 
ical temperatures, indicated by thin vertical lines, are ap- 
proached from above. In particular, we have T c ~ 0.43 T F for 
k F a = 1.9, T c ~ 0.36 T F for k F a = 1.8, and T c ~ 0.27 T F for 
k F a = 1.7. For k F a < ty/2 no ferromagnetism occurs within 
mean- field theory. In that case 1 / r s d (T) is smooth throughout 
the temperature range and exhibits the standard Fermi-liquid 
behavior l/r sd (T) oc T 2 for T -> 0. 



the enhancement of the conductivity in superconductors 
due to superconducting fluctuations above but close to 
the critical temperature [TT] . 

A closely related situation arises when the two layers 
of a 2D electron-electron bilayer are close enough to al- 
low the establishment of interlayer coherence [12] . In this 
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case, the two layers in the system can be labelled "up" 
and "down" along a "z"-axis, so that the which-layer de- 
gree of freedom becomes a spin one-half pseudospin. In- 
terlayer coherence in this language corresponds to pseu- 
dospin ferromagnetism with an easy x-y plane, since this 
orientation of the pseudospin describes a particle that is 
neither in the left nor in the right layer, but in a coherent 
superposition of the two. Furthermore, Coulomb drag 
becomes pseudospin drag, the mutual friction between 
two pseudospin states due to Coulomb scattering. This 
analogy prompted studies of spin drag, the frictional drag 
between electrons with opposite spin projection, in a sin- 
gle semiconductor [13] • While the realization of separate 
electric contacts to the two spin states remains an exper- 
imentally challenging problem, the spin drag is observed 
indirectly, by measuring different diffusion constants for 
charge and spin [14]. 

Because of the presence of other relaxation mecha- 
nisms, spin drag effects are usually not very large in 
semiconductors, and are even smaller in metals. This is 
completely different in cold atomic gases where scatter- 
ing between different hyperfine spins is the only mech- 
anism to relax spin currents, and was considered both 
for fermionic atoms [15] [16] , and for bosonic ones [17] . 
It is the purpose of this Letter to point out that a par- 
ticularly interesting situation occurs when spin drag is 
considered in a two-component Fermi gas that is close to 
a ferromagnetic instability [18-22 . Based on the anal- 
ogy between electron-hole bilayers and pseudospin ferro- 
magnets we expect that the spin drag will be strongly 
enhanced as the ferromagnetic state is approached from 
the normal side. Because atoms are neutral, the rele- 
vant experimental quantity is the spin drag relaxation 
rate, which for instance determines the damping rate of 
the spin dipole mode in trapped cold-atom systems [T6] 
and is thus accessible experimentally. Interestingly, an 
electronic analog of the spin dipole mode also exists [23] . 

Our main findings are illustrated in Fig. [I] This plot 
shows the spin drag relaxation rate l/r s d(T) as a func- 
tion of temperature, for various interaction strengths de- 
termined by the product of the Fermi wave vector k F 
and the scattering length a. The dramatic enhancement 
of the relaxation rate upon approaching the critical tem- 
perature for the ferromagnetic transition is clearly vis- 
ible. One of our motivations for considering this effect 
is the recent observation of ferromagnetic correlations in 
a two-component Fermi gas with strong repulsive inter- 
actions [24]. The fact that spin-polarized domains were 
not directly observed adds to the theoretical interest [25] 
in this experiment. The enhancement of the spin drag 
relaxation rate as the ferromagnetic phase is approached 
serves as a precursor probe for ferromagnetism that is 
distinct from, and adds to, the experimental methods of 
Jo et al. [24,, and is also interesting in its own right. In 
the following we present our calculations in detail, and 
present additional results and discussion. 



Spin drag relaxation rate. - We consider a 3D ho- 
mogeneous gas of fermionic atoms of mass m, with two 
hyperfine states denoted by | t) an d | |). The grand- 
canonical Hamiltonian with chemical potential \i is given 

by 



H = Jd 3 x J2 



2m 



U 



J d 3 x $\(x)$l(x)fa(x)fa(x) , (1) 



in terms of fermionic creation and annihilation operators 
^ a (x) and ifj a (x), respectively. At low temperatures s- 
wave scattering, described by U = 47raft 2 /m, dominates, 
and we have therefore omitted other interaction terms 
from this Hamiltonian. 

We first determine a frequency and momentum depen- 
dent scattering amplitude A^(q,uj) that takes into ac- 
count many-body effects on the scattering of atoms with 
opposite spin. We use the common random-phase ap- 
proximation that consists of summing all "bubble" dia- 
gram contribution to this effective interaction. This takes 
into account modifications of the interaction due to den- 
sity and spin fluctuations in an approximate way. The 
latter are essential as the spin fluctuations are strongly 
enhanced close to the ferromagnetic phase transition. In 
terms of the noninteracting (Lindhard) response function 
at nonzero temperature 



d 3 k 
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with e k = H 2 k 2 /2m and N k = [ e ^-^/ k ^ T + l]- 1 the 
Fermi-Dirac distribution function, the scattering ampli- 
tude reads 



U 2 3U 2 
= U + — Xnn(q^) — -r- 



xs z s z (q,u) , (3) 



where Xs z s z (nn)(Q,u) = Xo(<2,^)/[l ± Uxo(q,w)/2]- In 
this notation Xnn(Q^) is the density-density response 
function while Xs z s z (q^) describes the spin-spin re- 
sponse. The factor of three in the last term in the right- 
hand side of Eq. ([3| comes about because longitudinal 
and transverse spin fluctuations are both taken into ac- 
count. 

Within Stoner mean-field theory (MFT), ferromag- 
netism occurs when Xs z s z (0,0) diverges so that 1 + 
Z7xo(0,0)/2 = 0. This equation gives, together with 
the equation n = 2 J d 3 q N q /(27r) 3 for the total density 
determining the chemical potential, the critical tempera- 
ture T c as a function of k F a. For k^T c much smaller than 
the Fermi energy e F = h 2 k 2 ? /2m = ft 2 (37r 2 n) 2 / 3 /2m, this 
gives a critical temperature 
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FIG. 2: (Color online) Spin drag relaxation rate l/r s d as a 
function of kpa for various temperatures. The thin vertical 
lines indicate the critical values of the interaction parameter 
kpa at which ferromagnetism occurs. 



FIG. 3: (Color online) Panel a): spin diffusion constant D s as 
a function of temperature T for various interaction strengths 
kpa > 7r/2. Panel b): spin diffusion constant for k^a = 1.2 < 
tt/2. 



where v(sy) = m/cp/^ 2 ^ 2 is the density of states at the 
Fermi level. Note that one needs k^a > tt/2 for the 
critical temperature to be nonzero, and that there is a 
quantum critical point when k^a = tt/2 [T8] . 

Our next step is to use the scattering amplitude 
A^(q,oj) in Eq. ([3| in the well-known expression for 
the spin drag relaxation rate l/r s d, following from Boltz- 
mann theory p~3j [16] . This yields the result 
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sinh 2 [fiw/(2fc B T)] 



, (5) 



which, together with Eq. is the central result of this 
Letter and will be evaluated next. Before proceeding 
we note that our result is similar to the theory of spin 
diffusion in liquid He 3 [26], although the ferromagnetic 
phase transition was not considered in this context. 

Results. — In Figs, [I] and [2] we present the results of 
a numerical evaluation of Eq. ([5| , both as a function of 
temperature for various interaction strengths hp a, and 
as a function of interaction strength for various temper- 
atures. In experiments both dependencies can be ex- 
plored using a Feshbach resonance to tune the interac- 
tion strength [27 . Both figures clearly show the strong 
enhancement of the spin drag relaxation rate as the fer- 
romagnetic state is approached. The precise form of the 
enhancement is understood by keeping only the most di- 
vergent term in the scattering amplitude in Eq. ([3| so 
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with a(T) = 1 + f7xo(0,0)/2 - 1 - Uv(e ¥ ) + 
7r(T/Tp) 2 /12 + ... equal to zero at the phase transi- 
tion and positive for T larger than the critical tem- 
perature. In obtaining this approximation we have ex- 
panded around the ferromagnetic singularity and used 
that ^sm Xo(o,u) = — 7rv(sF)muj /2hk-Fq. Using these re- 
sults in Eq. (|5|, and expanding 1/ sinh 2 (x) ~ 1/x 2 , both 
the frequency integral and the momentum integral can 
be performed analytically if we use a cut-off of 2k^ on 
the momentum integration that diverges because of the 
expansion of 1/ sinh 2 (x). Ultimately we find in this man- 
ner that l/r sd (T) - l/rsd (T c ) oc (T — T c ) ln(T - T c ) for 
T I T c <C T F , which indeed accurately describes our nu- 
merical results near the critical temperature. 

We also consider the spin diffusion constant, which 
from the Einstein relation is given by D S (T) = 
cr s (T) /xs z s z (Q-> 0) an d the "spin conductivity" <r s (T) = 
nr s d(T)/m [28]. In Fig. [3] we show this constant as a 
function of temperature for various values of the inter- 
action strength. Near the critical temperature the spin 
diffusion constant vanishes as D B (T) oc (T — T C ) K with 
an exponent k = 1, because r s d(T c ) remains finite and 
Xs z s z (®i 0) diverges as 1/(T — T c ) within our random- 
phase approximation. 

At this point it is important to realize that from the 
point of view of critical dynamics our findings are mean- 
field like. If the spin dynamics can be effectively de- 
scribed by an isotropic Heisenberg ferromagnet (model J 
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[29]), the spin conductivity cr s (T) is expected to behave 
as £( 3-r 7)/ 2 very close to the transition, where the corre- 
lation length £(T) diverges as 1/{T — T c ) v and n and ^ 
are the usual static critical exponents of the ferromag- 
netic transition [29]. Since Xs z s z (Q,fy diverges as £ 2_?? , 
we find that the spin diffusion constant goes to zero with 
an exponent k = (1 — rj)v/2. In view of this possibility, 
we have therefore made sure, using the Ginzburg crite- 
rion [30] , that the upturn of the spin relaxation rate takes 
already place well outside the critical region where criti- 
cal fluctuations can be neglected and our random-phase 
approximation is appropriate. 

Discussion and conclusions — As we have mentioned 
in the introduction, the spin drag relaxation rate can 
be determined from the damping of the spin dipole 
mode [16] [23] in a trapped gas. Since the Fermi energy is 
usually much larger than the level splitting in the trap, 
our results show that the spin dipole mode is typically 
strongly overdamped, which makes the experiment more 
challenging. Nevertheless, such a measurement, as well 
as measurements of the spin diffusion constant as a func- 
tion of temperature, gives information on the proximity 
of the ferromagnetic phase transition. Although we have 
considered a homogeneous system, local-density approx- 
imations are generally valid for trapped Fermi systems, 
and in determining the damping of the spin dipole mode 
of a trapped two-component gas the homogeneous den- 
sity should in first approximation be taken as the central 
density of the atomic cloud. 

Within our present approach, we consider the transi- 
tion to ferromagnetism within MFT, which predicts it 
to be continuous. One interesting aspect is that taking 
into account correlation effects beyond MFT [22] results 
in: (i) an increase in the critical temperature for a given 
value of k^a and (ii) a change in the character of the 
transition from second to first order at very low tem- 
peratures. We expect that these modifications will not 
qualitatively affect the upturn of l/r s d, as long as one 
remains outside the critical-fluctuation region. 

In future work we intend to explore the effects of lower 
dimensionality, and the implications of critical and quan- 
tum critical fluctuations on the exponent k that deter- 
mines the behavior of the spin diffusion constant. Fur- 
thermore, as our approach is distinct from the work of 
Hu [10] , who considered the upturn of the Coulomb drag 
resistivity as one approaches the exciton-condensed state 
in bilayers, we intend to study this system as well with 
our present approach. 
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